A quantum phase transition in a quantum external field: Superposing two magnetic 

phases 
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We study an Ising chain undergoing a quantum phase transition in a quantum magnetic field. Such 
a field can be emulated by coupling the chain to a central spin initially in a superposition state. 
We show that - by adiabatically driving such a system - one can prepare a quantum superposition 
of any two ground states of the Ising chain. In particular, one can end up with the Ising chain in 
a superposition of ferromagnetic and paramagnetic phases - a scenario with no analogue in prior 
studies of quantum phase transitions. Remarkably, the resulting magnetization of the chain encodes 
the position of the critical point and universal critical exponents, as well as the ground state fidelity. 

PACS numbers: 



Quantum phase transitions (QPTs) occur when dra- 
matic changes in the ground state properties of a quan- 
tum system are induced by a tiny variation of an external 
parameter, such as the magnetic field in spin systems [l| 
or the intensity of a laser beam in cold atom simulators of 
Hubbard-like models [2]. In all current studies of QPTs, 
the external parameter is assumed to be classical, i.e., 
it has a well-defined instantaneous value. However, the 
field inducing a QPT can be quantum as well, taking 
on different values by virtue of being in a superposition 
of states. In fact, tremendous progress with the prepa- 
ration and manipulation of cold atom/ion systems will 
allow for creation of scenarios where the quantum nature 
of the "external" parameter will play a central role. 

For instance, cavity- QED systems offer intriguing pos- 
sibilities to study quantum control parameters [3|-[5|. In 
these systems, photons - bouncing off two parallel mir- 
rors - interact with ultracold atoms. If the number of 
photons in the cavity does not fluctuate, atoms expe- 
rience an "external" periodic potential cos^(kx), whose 
amplitude is proportional to the number of intra-cavity 
photons (k is the photon wave- vector). Atoms in such a 
system would be either in the superfluid phase or in the 
Mott insulator phase [2]. It may be possible, however, to 
create a coherent superposition of the intra-cavity pho- 
tonic states, giving rise to quantum fluctuations in the 
number of photons between the mirrors. The atoms 
would then be exposed to a coherent superposition of 
periodic potentials with the same period but differing 
amplitudes. In this case, one can have atoms in a su- 
perposition of two quantum phases, i.e., simultaneously 
in superfluid and Mott insulator ground states [3]. Such 
a situation has no counterpart in traditional studies of 
QPTs where the system is either in one phase or another. 

An analogous phenomenon can be envisioned in central 
spin models. These models are used to describe qubit 
- environment interactions in nitrogen-vacancy centers 
in diamond [6], quantum dots in semiconductors [7, 8], 
NMR experiments j^], etc. The focus is typically on the 
loss of coherence of the qubit while ignoring the environ- 
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FIG. 1: Schematic of a superposition of two different quan- 
tum phases in an Ising chain (in the paramagnetic phase spins 
try to align with the classical magnetic field, which is oriented 
in the z direction here; in the ferromagnetic phase spin inter- 
actions try to align spins perpendicular to the field; see the 
Hamiltonian in Eq. ([T]) for details). One can prepare such a 
state by adiabatically evolving the chain in the presence of a 
central spin followed by a measurement of the spin. 



mental degrees of freedom. We will take the opposite 
perspective and explore the quantum state of the en- 
vironment subjected to an effective quantum potential 
originating from the central spin. For an experimental 
study of such a scenario, one needs a well-controlled sys- 
tem, which we expect will be delivered in the foreseeable 
future by ion simulators of spin chains [10|. 



I. RESULTS 

The model. We will discuss the most striking conse- 
quence of a QPT in a quantum potential: The possibility 
of having the system in a superposition of ground states 
belonging to different phases, as shown in Fig. [TJ We 
consider a quantum Ising chain uniformly coupled to a 
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FIG. 2: (a) Schematic of the central spin model in a classical 
magnetic field: the central spin is equally coupled to all the 
spins- 1/2 arranged on a ring, (b) Possible realization of the 
central spin model in a linear ion chain. The ions emulate 
the effective spins- 1/2. The couplings between the effective 
spins- 1/2 are optically engineered to be the same as in the 
(a) panel. One of the ions is differently coupled to the rest of 
the chain to play the role of the effective central spin. The 
effective magnetic field is also optically engineered. 



(central) spin- 1/2 (Fig. EJ: 
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where ^ 1 is the number of spins arranged on a peri- 
odic ring. The central spin contribution is contained in 
the effective magnetic field operator 
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where g is the (classical) magnetic field strength and Sa^ 
is the quantum component of the field generated by an 
Ising coupling to the central spin (0 < <C 1). Without 
the coupling to the central spin, i.e., when H{g, S = 0) = 
Hi{g), the Ising chain in the ground state is either in 
the ferromagnetic phase (|^| < 1) or in the paramagnetic 
phase (1^1 > 1), with critical points at = if- 

The successful implementation of a recent proposal 
simulating arbitrarily-connected spin models in linear ion 
chains [10] will put the system we consider within exper- 
imental reach. It is interesting to note that there is no 
need to arrange ions, i.e., effective spins, on the ring to 
simulate the Hamiltonian ([1]), see Fig. [2fa). So far the 
simulation of an A/" = 6 Ising chain with long-range inter- 
actions between the effective spins- 1/2 has been demon- 
strated [11], analogous to that shown in Fig. [2fb). The 
proposal is scalable and it is expected to allow for quan- 
tum simulation of models with N ^ 1 effective spins. 

Preparation of the superposition state. QPTs 
can be studied either by diagonalizing the Hamiltonian 
for a fixed set of coupling parameters or by adiabati- 
cally evolving the system from an easy-to-prepare ground 
state, which is especially relevant in cold atom experi- 
ments [2,]. We take the latter approach, as the former 



will always force the central spin to point in either -\-z or 
—z direction for any g ^ because [H, erg] = 0. 

We assume that at t = the chain is prepared in a 
ground state and its coupling to the central spin is turned 
off, which provides freedom to engineer the state of the 
central spin. The composite wave function is \tlj{g{ti))) = 
\S)\g{ti)), where \g) is a ground state of Hi{g) and the 
central spin state is 



|5)=cte'^t(«.)||)+c;e 



where Cf^ > 0. By changing both the bias field g and 
the coupling 5, the wave-function evolves according to 



mg{t)))=fexp 



i-i: 



dtH[g{t),5{t)]]mi)) 



where T is the time-ordering operator. 

As was shown in Ref. [12], \tlj{g{t))) can be simplified. 
Considering adiabatic evolution, we obtain 

= e'^^^'^ct It) \g + S)+ e^^^(*)c; \i) \g - S). (3) 

We thus study finite, i.e., gapped, systems so that adia- 
batic evolution is possible by changing g{t) and S{t) slow 
enough. In the state (|3j), the chain experiences an aver- 
age magnetic field 
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with fluctuations 



In particular, this shows that once the desired coupling 
5 is adiabatically reached, fluctuations of the quantum 
potential are fixed. 

The state ([3]) is already a Schrodinger's cat state, 
where the two "macroscopically" distinct possibilities 
are the ferromagnetic and paramagnetic phases, both 
of which are coupled to the auxiliary two-level system. 
Since we are interested in the QPT of the Ising chain, we 
will "trace out" the central spin by measuring its state. 
If we will do the measurement in the , \\)} basis, the 
superposition will be destroyed and the state of the Ising 
chain will be one of the ground states \g ^5). Measure- 
ment in any other basis will result in a superposition of 
Ising ground states at different magnetic fields. 

We assume that the measurement will be done in the 
eigenbasis of the operator: 
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In this basis, 
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where we write (j)^^^ as shorthand for ^^^^ (t). There- 
fore, the measurement of the central spin in the state |±) 
leaves the chain in the state 



^\ ± 2ctc^ cos (0t - 04.) T{g, S) ' 



(4) 



where T{g^ (5) = — + is a ground state fidelity |13j , 
or simply fidelity, whose crucial role in this problem will 
be carefully discussed below. Without loss of generality, 
we define it in such a way that F{g^5) > 0. 

Now we comment on the measurement of the central 
spin. The state |±) will occur with probability 

P± ih - ^i) = ^ ± c^c; COS {(j)^ - 7", 

which depends on the relative phase between the ground 
states in the superposition (|3|). Since the point of the 
measurement is to prepare a well-defined superposition 
state of the Ising chain, we will describe properties of the 
Ising chain after finding the central spin in, e.g., the |+) 
state. Then, the Ising chain will be in the state 



I Ising) 



V'l + 2c^ci cos ((/)^ - J^{g, 8) 



(5) 



This state is the desired superposition of ferromagnetic 
and paramagnetic ground states when 

g-8<gc<g^8, 

which is depicted in Fig. [TJ We propose to call such 
a state a Schrodinger magnet. The possibility to create 
such a novel state of matter is offered by the quantum 
magnetic field in Eq. (j2j). Indeed, if there would be no 
quantum component in the magnetic field, the wave func- 
tion of the Ising chain after the adiabatic evolution would 
correspond to either a ferromagnetic or a paramagnetic 
phase ground state, but never to a superposition of both. 

Ising chain in the superposition state. For sim- 
plicity, we assume that the measurements on the Ising 
chain are performed immediately after measuring the 
central spin. The expectation value of an operator O 
in the state (|5]) is 



O = (lsing|0|lsing) = 



+ 2c^c; cos (A) 0+- 
1 + 2c^ci cos (A) T 



(6) 



where A = (p-^ 



and 



ciO++ 



c^O- 



{g±5\d\g±5), 



is the "standard" average, and 

0+- = {g + S\d\g-S) 

designate the cross term that arises. For clarity of 
presentation, we restrict ourselves to real 0+~ in Eq. 



(j6j), because 0+ is always real for the operators O = 
cr^, cr^, cr^cr^+i that we study. It is not real, however, 
for all operators (e.g., for O = a^), and it is a straightfor- 
ward exercise to extend our calculations to these cases. 
While the standard average does not yield any new infor- 
mation, the cross term provides a non-trivial correction 
absent in a quantum phase transition in a classical field. 

To further simplify the discussion, we average Eq. ([6]) 
over several realizations where the appearance of the rel- 
ative phase A of the superposition (|3|) is given by some 
probability distribution p{A). For example, such aver- 
aging may appear due to preparation of the central spin 
with random initial phases <p^^i (ti). We assume below for 
simplicity that p (A) = ^ for A G [0,27r). We denote 
the result of such averaging as O, and define its variance 



through var(O) = 
notation 



O . Finally, we introduce the 



0+- = 0+-/^, 



because for the O operators that we consider, is 
a well-defined non-zero quantity in the thermodynamic 
limit in which fidelity typically tends to zero (see the 
Discussion section). 

The phase- averaged observable and its variance are 



-_ /,-^"dAp(A)P+(A)0 
/o'^dAp(A)P+(A) 
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with X = 2c^c^T . By expanding 
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we see that the square root of variance is proportional to 
fidelity when x <C 1. The role of fidelity in our problem is 
discussed in the Discussion section. In the following, we 
use the exact solution of the Ising model to study expec- 
tation values of different observables in the superposition 
state ([5]), see the Methods section for details. 

We start by looking at O = = cf^. Mf^ terms 
have been calculated in Ref. 
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where x± = ^{d ^ ^)/(l + ^ ^ i ^^^^ K and E are 
elliptic functions of the first and the second kind, respec- 
tively. Above a large limit is assumed to simplify the 
expressions (see the Methods section for the exact finite 
N expressions). 

The cross terms can be obtained from the eigenequa- 
tion 



Hi{g±5)\g±5)=Ne{g±5)\g±5), 



(7) 
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FIG. 3: Mean magnetization along the z-direction in a quan- 
tum superposition of different phases, O = Mz = cr^. The 
plot shows exact results obtained from expressions listed in 
the Methods section. The solid line shows Mz evaluated from 
Eq. (|6|), the dashed line shows Mz = M|, and the dashed- 
dotted line shows Mz =1= ^var (Mz). We assumed N = 100, 
6 = 0.05, = 1/2, = -^3/4, and g — 1. For these parame- 
ters J" ^ 0.41, Ml 6.0 X 10"\ and ^var(M^) 8.9 x 10"^ 



where e is the ground state energy per site. Indeed, one 
gets from it 



M+- =T{g,S) 



s{g -S)- e{g + S) 
25 



In the limit of AT ^ oo, £(5 ± <5) = + 5 ± 5\E{x±) 

[T3 |. Consequently, 
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The dependence of magnetization on the relative phase 
of the superposition and the variance of magnetization 
at = 1 are depicted in Fig. [H 

We mention in passing that similar expressions can be 
obtained for O = Cx = cr^cr^+i- Indeed, it is known from 
Ref. y^] that 
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and one can use again Eq. ([7]) to derive 



K{x±), 



c+- = 



^-^e{g + 5) - ^-—e{g - 5) 



Since these results are analogous in structure to the ones 
already discussed, will not analyze them. 

Next, we study spontaneous magnetization in the x- 
direction. The system will acquire such a magnetization 
when a tiny field breaking the —a^ symmetry of 

the Hamiltonian is present. When necessary, we thus add 
a —hY^^^^ (7^ term to Hi{g) and denote a ground state 
of the resulting Hamiltonian as \g, h). Without the quan- 
tum magnetic field, ^ = 0, the Ising chain acquires macro- 
scopic magnetization (along the direction of the symme- 
try breaking field h) only in the ferromagnetic phase. 



FIG. 4: Spontaneous magnetization in the x-direction in a 
quantum superposition of different phases. The red dashed 
line is the "standard" average = 4 M++ + cl M--. The 
shaded area between the solid green lines marks the range of 
variation of Mx due to variation of the relative phase 0^ — 0^ 
(similar variation, but at a single magnetic field is de- 
picted in Fig. [3]). This is a numerical result obtained with 
the "periodic" TEBD algorithm for 6 = 0.05, h = 0.0001, 
= = 1/v^, N — 100, and x = 50 (the cut-off param- 
eter of the algorithm). The spontaneous magnetization does 
not disappear for ^ > 1 + 5, when both states in the super- 
position are in the paramagnetic phase, due to the non-zero 
symmetry-breaking field h (see Fig. [5] for the h ^ 0^ and 
1 limits; see the Methods section for details). 



This magnetization can also be calculated by studying 
the correlation function jl4| 
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'Wf(^R\9)\ = \im {g,h\cT^\g,h) = (1 



2)1/ 



Importantly, it encodes the critical exponent /3 = 1/8 
115|. 

To study spontaneous magnetization in the presence 
of the superposition of ground states, we find numer- 
ically the states \g ± S^h) using a periodic version 
of the TEBD algorithm [13, IHj. Then, we cal- 
culate M±± = {g±S,h\a^\g±6,h)\^^Q and M+- = 
{g -\- S^h\a^\g — 5^h)\^^Q. Naturally, for large enough 
systems, the standard result is reproduced by numerics: 



±± 



211/ 



for 1^ ± (5| < 1 and zero otherwise. The results of TEBD 
calculations are plotted in Fig. |4j The presence of the 
cross term magnetization, resulting from the superposi- 
tion of two ground states in Eq. ([5]), leads to sizable 
deviations from the "standard" average. 

To analyze this deviation more efficiently in the ther- 
modynamic limit, we study the asymptotic behavior of 
the two-point correlation functions: 



lim 

R^oc 



l+R/T 



(8) 



where (• • • = {g-\-S\ • • • \g — S)/J^. It can be done using 
the exact solution of the Ising model through fermioniza- 
tion, where we express the correlator as a determinant of 
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FIG. 5: Scaling properties of the spontaneous magnetization 
cross term in the x-direction (|9|). Upper panel: Illustration 
that M^-^ {g, 6) at the critical point g — Qc — 1 scales as 5^^^ . 
Crosses show numerics based on Eq. (|8|), while the straight 
line is a fit. The same result is obtained near the critical point 
for other c = [q — Qc) and ^ < 1 + ^. Lower panel: Illustra- 
tion of the scaling function B{c) near the critical point and 
far away from it. Crosses show numerics, the solid line con- 
nects them, and the dashed line is (— 2c)^/^. See the Methods 
section for details. 

a 2R X 2R block Toeplitz matrix, which is then numeri- 
cally evaluated (see the Methods section for details). 

As shown in Fig. O we find that the scaling of 
around the critical point is consistent with the ansatz 

Mt-=F5^B{c), c=^, (9) 

where /3 = 1/8 and Qc = 1 for the Ising chain that we 
study, and B{c) is the scaling function. It is nonzero 
when at least one of the superposed ground states is in 
the ferromagnetic phase, i.e., B{c) ^ for c < 1. Far 
away from the critical point, we observe that ^ 
^ M^- and so B{c <C -1) (-2c)i/^ 

II. DISCUSSION 

We have seen that the presence of the quantum ex- 
ternal field allows for creation of the superposition state 
of two distinct ground states in general and two distinct 
phases in particular. If this happens, expectation values 
are altered by the cross terms. The magnitude of this 
effect can be sizable as depicted in Figs. [3] and HI 

A fundamentally important question can now be an- 
swered: what is the role of the system size in a quantum 
phase transition in a quantum field and what critical in- 
formation is imprinted onto the cross terms. 

To answer it, we note that all the cross terms that we 
studied are a product of the two terms: the ground state 



fidelity F{g^5) and a term that has a well-defined non- 
zero value in the thermodynamic limit. Ground state 
fidelity, however, typically disappears in the thermody- 
namic limit of A/" ^ oo often invoked in the context of 
quantum phase transitions. 

This is known as Anderson orthogonality catastrophe 
after the seminal work reported in Ref. [19]. Therefore, 
we are interested in the studies of systems for which N ^ 
1 (to see quantum criticality ) , but still A/" < oo (to avoid 
the catastrophe). There are three options here, which we 
will discuss below. Instead of providing specific results 
for fidelity of the Ising chain, we provide general scaling 
results to highlight the role of critical exponents in our 
problem and to keep the discussion concise. 

First, one can consider the limit of S ^ taken at 
fixed AT > 1. Then, fidelity reads [20H22| 

\nJ^{g,5) - (10) 

near the critical point. Here, d is system dimensionality 
and u is the critical exponent (correlation length diverges 
as \g — gc\~^ near the critical point; d = u = 1 in the Ising 
chain that we consider). Since fidelity is close to unity 
in this limit, the cross terms do not get small. One must 
remember, however, that they will be dominated by finite 
system size corrections requiring a separate study, which 
is beyond the scope of this work. 

Second, in the limit of A" ^ oo at fixed S - the one 
that we assumed in our calculations - one can focus on 
the "moderately" large systems. To explain this term we 
note that near the critical point (in the above-mentioned 
limit) m 

\nJ^{g,S) - -NS'^^ (11) 

The crossover from Eq. ([TQ|) to Eq. (pT]) happens near 
the critical point when |23| 

A^^^^ - 1. 

We define the "moderately" large system to be just large 
enough to exhibit the scaling of fidelity with S and A^ 
given by Eq. (fTTj) rather than Eq. (fTQ|) . In the Ising case, 
Eq. (dl]) predicts InT -NS while Eq. ([10|) predicts 
In ~ —N'^S'^. For a "moderately" large system fidelity 
shall not be too small to erase the contribution of the 
cross term (see, e.g.. Fig. H]). 

Third, one can study superpositions of two ground 
states from the same phase far away from the critical 
point. There In J* - -A"^^!^ - ^c|'^''"^ and for A" > 1 
fidelity can be kept close to unity by a proper choice of 
6. The downside of this scenario is that we loose the 
possibility to superimpose two phases. 

From the above discussion, we see that the critical ex- 
ponent u is imprinted onto the cross term via fidelity. 
Also the critical exponent P is seen in the cross term 
contributing to spontaneous magnetization. The 
location of the critical point is most directly seen in the 
cross terms M+~ and through "divergence" of their 
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second derivative over g taken at Qc ± 5. This is caused 
by the singularity of the second derivative of the ground 
state energy per spin across the critical point. This sin- 
gularity will be rounded off in finite systems {N < oo), 
but nevertheless there shall be pronounced peaks visible. 
We also note that while the "standard" averages M| and 
(7| also encode the position of the critical point, they do 
not encode the critical exponent u. 

To observe the superposed phases, experiments will 
have to keep decoherence to a minimum. The decoher- 
ence rate of the state in Eq. (jSj will depend on how 
well the environment distinguishes the two components, 
which will depend on the system size (see, e.g., Ref. [24]) 
and the overlap between the two states (fidelity). Thus, 
N cannot be too large and T cannot be too small. This 
is a similar issue to being able to observe the effect of 
the cross term, which we discussed above. We thus do 
not expect decoherence to be overwhelming in a properly 
prepared setup. Further, the system size can be used as 
a parameter controlling decoherence, and its manipula- 
tion should be sufficient to bring decoherence down to 
an acceptable level. Looking from a different perspec- 
tive, studies of decoherence of such a novel macroscopic 
quantum superposition are fundamentally interesting on 
its own, e.g., to boost understanding of the quantum-to- 
classical transition. 

To conclude, we considered a quantum phase transi- 
tion of an Ising chain exposed to a quantum external 
field. This scenario can be used to create a new state of 
matter where the system is simultaneously in two distinct 
quantum phases. Observables on the chain then take on 
forms that encode the ground state fidelity, the location 
of the critical point, and the universal critical exponents 
of the system. These findings set the foundations for de- 
veloping a scaling theory of quantum phase transitions 
in quantum fields. Recent advances in cold atom cavity- 
QED and ion traps may lead to experimental realization 
of superposed phases. 

This work is supported by U.S. Department of Energy 
through the LANL/LDRD Program. MMR acknowl- 
edges support from the FWF SFB grant F4014. 



III. METHODS 



Here we provide some technical details regarding our 
calculations. 

The Ising Hamiltonian Hi{g) is diagonalized using the 
standard approach (see e.g. Ref. [25]). The Jordan- 
Wigner transformation, 

< = 1 - 2c];Cn, < = (Cn + 4) n ~ 24c^), 



where Cn are fermionic annihilation operators, transforms 
the Ising chain to a free-fermion model. After applying 



the Fourier transform 



-i7r/4 



the Hamiltonian takes the form: 

Hi{g) = ^(2^4 - 1) (^ - cos k) + (4cL/e + C-kCk) sin k, 

k 

IT N 
/C = ±(2.+ 1)-, , = 0,...,y-l. 

Diagonalization of the Hamiltonian with the help of the 
Bogolubov transformation leads to the following ground 
state wave-function 

1^ ± ^) = n [^^^ i^kM \^kO-k) - sin (0^/2) \lkl-k)] , 

k>0 

where \mk^m-k) describes the state with m = 0, 1 pairs 
of Ck quasiparticles with momentum k and 



tan 0^ 



sin k 



g ± S — cos k ' 



To prepare Fig. [3l we fix the system size N and use the 
following exact expressions for magnetization and fidelity 



{5±<5|<|5±<5> = ^^cos0± 

k 

COS 

N 



J" 



-S\g-5) = H' 



k COS ^ 2 ^ 



k>0 



> 0. 



To prepare Fig. IH we select ^, S and /i, and calculate 
the ground states \g ± d^h) of the Ising chain exposed 
to transverse and longitudinal magnetic fields. This is 
done through imaginary time evolution performed with 
the periodic TEBD algorithm. A global phase of the 
wave- functions is then chosen to make = {g — S^h\g 
J, h) positive. We then directly calculate and 
(both are positive). Putting these results into Eq. (j6j), 
one can calculate the spontaneous magnetization in the 
x-direction in the superposition state (j5j). The result 
is still dependent on the relative phase (j)^ — (j)^. When 
this phase is either or ±7r, spontaneous magnetization 
at any fixed g, ^, and h reaches an extremum. These 
extremal values are depicted by solid green lines in Fig. 



m 



To prepare Fig. O we calculate correlation function 



where we introduce bn = c}^ — Cn and 



Cn. We 



study it, because Mj^ = limi^^oo \jctx{R)T' 

The next step is to use Wick's theorem extended to 
such a cross-correlation [26]. It can be used as long as 
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the overlap T is nonzero, which is the case in our cal- 
culations. Then extending the results of Ref. [27], we 
find that C^~{R)jr can be expressed as a Pfaffian of a 
2R X 2R antisymmetric matrix, which can be converted 
into a determinant: 



C+-(i?)^= x/det [Ar], 

_(^m+l^n)J (^m+l^n+l)J 



Ar 



m,n=l...R 



where Ar is a block Toeplitz matrix. Apart from a few 
special cases, it is not known how to calculate such a de- 
terminant analytically [28]. Thus, we use numerics with 
a large enough R to obtain a well-converged result. We 
employ a continuous (i.e. N oo) approximation for 
the elements of the Toeplitz matrix 



- r 



dk tan^ — ^ gifc(m-n) 



-1 
2^ 



dk 



-iie++e-)/2 



^ik{ra — n) 



and {amCLn)^ = (bmbn)^ , {amK)jr = -{Kam)jr • 

Regarding the parameter i?, we mention that it has to 
be of the order of 500 (2000) for g = 0.995 and S = 0.01 
{g = 1.005 and S = 0.01) in order for the results to 
be converged to the R ^ oo limit. For every g and 5 
sufficiently large R is chosen to calculate data for Fig. [5l 

We also mention that we verified the Pfaffian-based 
numerics with a direct numerical calculation using the 
TEBD algorithm. For systems composed of about 100 
spins, for which the TEBD algorithm can still be ef- 
ficiently applied, spontaneous magnetization from both 
calculations agree. 

Finally, we provide definition of the elliptic functions 
that we use in the Results section: 



.7r/2 

K{x) = / — 
^0 VI 



X sm 



^7r/2 I 

E{x) = / d(j)\/ 1 — X sin^ ^. 

^0 



cos ■ 
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